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t^J- , to investigate the importance of final state interactions [FSI] at large momen- 

turn transfers in inclusive quasielastic electronuclear reactions. The optical 
potential is derived from first-order multiple scattering theory and then is 
used to calculate the FSI in a nonspectral Green's function doorway approach. 
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Abstract 



A relativistic finite nucleus calculation using a Dirac optical potential is used 



At intermediate momentum transfers excellent predictions of the quasielastic 



56 Fe(e, e') experimental data for the longitudinal response function are ob- 
tained. In comparisons with recent measurements at \q\ = 1.14 GeV/c the 
theoretical calculations of Rl give good agreement for the quasielastic peak 
shape and amplitude, but place the position of the peak at an energy transfer 
of about 40 MeV higher than the data. 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



There has been a recent interest in inclusive quasielastic electronuclear reactions at large 
momentum transfers (|g|>l GeV/c), especially with the upcoming programs planned at 
CEBAF. Along with recent measurements [l|] there has been some theoretical work investi- 
gating the physical role of relativity and final state interactions [FSI] in such reactions . 
From these studies it has been noted that a full consistent finite nucleus calculation would 
be helpful in discerning the various physical contributions. This paper attempts to address 
that need. 

In the simple relativistic Fermi gas calculation of Ref. |2| the implications were that the 
role of FSI at |g|> 1 GeV/c appeared to be greatly reduced, especially since the use of a 
real energy-independent potential to model the FSI caused the predicted position of the 
quasielastic peak to move significantly away from the data. The use of a Fermi gas model 
may be misleading however, since the recoil effects will be misrepresented which may affect 
the calculated position of the quasielastic peak. A finite nucleus model should be more 
appropriate. In Ref. an energy-dependent real potential was introduced to give the FSI, 
and the conclusion from this study is that the FSI remain important at large \q\, and that the 
energy dependence is required to enable one to predict the peak position correctly. In that 
work the imaginary part of the optical potential is neglected. From optical model studies 
of elastic nucleon-nucleus scattering and from multiple scattering theory there is known to 
be a strong energy dependence in the optical potentials. To delete the imaginary part of an 
energy-dependent potential is to break the unitarity constraint and thus incorrectly represent 
the reactive content of the optical potential. For a more physical representation, the full 
complex energy- dependent optical potential should be included in a consistent manner. The 
importance of such considerations was discussed in detail in Ref. |4j] . 

In this paper two complex energy- dependent Dirac optical potentials derived from mul- 
tiple scattering are used in a relativistic finite-nucleus calculation to calculate the separated 
response functions for inclusive quasielastic (e, e') scattering from 56 Fe. In these calculations 
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Dirac dynamical effects resulting from couplings to negative energy states, which were shown 
in Ref. Q to be important, are included. In Section II a theoretical discussion of the model 
and the calculation is presented. The results and the comparison with the experimental 
measurements are presented in Section III followed by a conclusion. 



II. THEORETICAL DISCUSSION 

Although the main mechanism in quasielastic reactions is assumed to be the knockout 
of a single nucleon, in inclusive reactions all possible final states are included in the ex- 
perimental measurements. The optical potential implies the existence of other final states 
besides the knockout channel within the imaginary part. For this reason to include all of 
the possible final states implied by an optical potential the nonspectral Green's function 
doorway approach ||, which is discussed more fully in Ref. ||, is used. The longitudinal 
and transverse response functions within the one-photon-exchange approximation are given 
by: 

R L (q,cu) = W 00 (q,u) 

R T (q,co) = W u (q,u) + W 22 (q,u;), (2.1) 

where 

= T$ (ilWl/W + w - E f )(f\r(q)\i). (2.2) 

i f 

Here \i) represents the initial nuclear many-body state, while the sum over |/) corresponds 
to all final states of the full hadronic many-body assembly. J^{q) is the electromagnetic 
nuclear current operator and the J2i denotes an average over the initial states. 

One would like to perform an explicit sum over the complete set of complex inelastic 
reaction channels in the final state, but in practice such a many-body calculation is pro- 
hibitively difficult. Therefore a nonspectral approach is used where the sum over all of the 
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final states within a particular space is implicitly performed, by considering the full A-body 
Green's function. 

Suppressing the discrete state contribution, eq. fl2.2| ) can be rewritten in terms of the 
forward virtual Compton amplitude: 

W^(q, u) = --ImT^iq, to), (2.3) 

7T 

where 

T^(q,u) = ^\n q yG(uj + E^riq)^). (2.4) 

i 

Here G is the full many-body propagator for the A— nucleon system. If J M is assumed to 
be a one-body operator, it can be shown that within this one-body space G reduces 
to the optical model Green's function. By using the optical model Green's function in a 
nonspectral form, then a proper and consistent unitary description of the reactive content 
of this inclusive reaction is maintained 0]. If one assumes that each knockout channel is 
represented by the same optical model potential, then the following substitution can be 
made: 

G ~ G op t — Go + Go U op t G op t- (2-5) 

G op t corresponds to the use of an optical model potential to represent the final state inter- 
actions between the ejected nucleon and the residual nucleus. G is the free propagator for 
a nucleon within the nuclear medium. 

To reduce the calculation to the plane wave approximation [PWA], U opt is set to zero or 
equivalently: 

G « G , (2.6) 

which leads to eq. (|2.2|), where only final plane wave states for the ejected nucleon are 
considered. 

The calculations are performed in a fully-off-shell momentum space representation: 



T^(<7» = E/ - j y L (Ap-$n-q)G opt {pJ,E)r{p')(p -q\t) . (2.7) 

i 

Here G opt (p, if; E) is calculated from the fully off-shell relativistic optical potential as derived 
from multiple scattering theory. The optical model Green's function is calculated as the 
solution of the Lippmann-Schwinger equation in momentum space to give the fully off-shell 
nucleon-nucleus T-matrix. The equations are solved in partial wave form so as to include 
the spin-orbit contributions in a convenient way. 

G opt (pJ; E) = G (p; E) 6®(p- f) + G (p; E) V opt (p,f) G opt (p,if; E) (2.8) 

= G (p; E) 5^(p- p) + G (p; E) T opt (p,jf; E) G (p'; E) 
T opt (p, if; E) = V opt {p, p) + J d 3 f V opt (p, f) G (p"; E) T opt (f, f; E) (2.9) 

In the above equation the free Dirac propagator can be separated into a positive-energy 
projecting part and a negative-energy projecting part, so that contributions to eq. Q2.7D 
involve coupling to the negative-energy Dirac sea. It has been shown that such effects can 
play a major role in the calculated response functions 

In the PWA and FSI calculations a single particle description is used. Bound state 
wave functions are taken from a Dirac-Hartree calculation and are represented in Dirac 
4-spinor form. Since the Dirac-Hartree calculation performed in Ref. |J assumes spherical 
symmetry and 56 Fe is not a doubly magic closed shell nuclei, approximations are used to 
represent the valence nucleons. In this case the valence shell is represented as a closed shell, 
but with fractional occupation numbers to give 26 protons and 30 neutrons. The current 
operators are treated in relativistic form and no nonrelativistic reduction is performed. 

The analysis presented in this paper will be performed for both a.) the relativistic plane 
wave approximation and b.) with the Green's function doorway approach to include the FSI. 
The form of the electromagnetic nucleon current inside of the nucleus is unknown, hence 
two different functional forms of the free electromagnetic current operator are used: 
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and (2.10) 
2m 

where G m = Fi + F 2 is the familiar Sachs magnetic form factor and K, = k + k'. The bars over 



JC and q indicate that k and k' are fixed to the onshell values, e.g. k° = ±y k 2 + m 2 , where 
the sign is dependent upon the (±) energy character of the Dirac spinor. The definitions of 



J CC 2 and J cc i correspond to the cc2 and ccl operators defined in Ref. ||10|| . These two operators 
are identical on-shell and hence give the same free nucleon electromagnetic representations. 
In the off-shell case where one scatters from bound nucleons or when one includes FSI, these 
operators give differing results. The most general form of the current operator contains 12 
independent terms, in which only two independent terms survive in the on-shell limit. To be 
able to construct the complete operator with the accompanying form factors would require 
a reliable off-shell nucleon structure model, for example a QCD based model. Hence the 
calculated differences between the ccl and cc2 operators can only be understood in terms 
of the underlying nucleon structure. For a detailed analysis of the effects and uncertainties 
represented by these two operators, please see Ref. [f^ |lllJl2"| . 

Current conservation is imposed by means of the standard replacement of q- J by q J°/\q\- 
In general the current is not conserved by these two operators, since there typically is not a 
consistent Hamiltonian treatment of the initial state, the final state and the electromagnetic 
current interactions. The form factors used in this paper are taken from Ref. |13| . 



Two relativistic complex optical potentials derived from multiple scattering theory are 
used to represent the FSI. In this case nucleon-nucleon [NN] t-matrices are folded with local 
densities in the optimum factorization approximation [14j] to give the optical potential. The 
negative-energy part of the optical potential is constructed using an approximate approach 



T5| . For one optical potential the NN t-matrices are calculated from the fully-off-shell 
full Bonn potential fl6| , which include the effects of relativistic kinematics, retarded meson 
propagators as given by time-ordered perturbation theory, and crossed and iterative meson 



exchanges with NN, NA and AA intermediate states. For ejectile energies greater than 
300 MeV, an extension of the Bonn meson exchange interaction above pion production 
threshold is used |T7| . The second optical potential uses the NN interaction of Ref . |jTB| . 



The proton densities are taken from electron scattering measurements [OJ, while the neutron 
densities are those calculated from the Hartree-Fock-Bogolyubov calculation of Ref. [20]. In 
calculations of elastic nucleon-nucleus scattering, the use of the Bonn potential tends to give 
a better representation of the data than the Franey-Love amplitudes, probably due to the 
superior off-shell behavior of the Bonn potential. 

In this paper the nonrelativistic calculations of FSI are constructed from the relativistic 
calculations with the exception that all of the negative energy contributions which result 
from the Dirac dynamics are neglected. This includes those negative energy contributions 
that arise from the construction of U opt and in the calculation of G opt . This is the manner 
in which the nonrelativistic calculation is calculated here, where relativistic kinematics are 
maintained. 



III. RESULTS AND COMPARISONS WITH THE DATA 

To gauge the accuracy of the theoretical model, comparisons with quasielastic 56 Fe(e, e') 
data are made at \q\ = 410 and 550 MeV/c. These comparisons are made both with 
nonrelativistic FSI and with relativistic FSI including Dirac dynamical degrees of freedom. 
The PWA results with no FSI (V op t = 0) and the nonrelativistic FSI calculations are shown 
in Figs. 1 and 2 for \q\ = 410 and 550 MeV/c, respectively. The PWA calculation places the 
peak position for Ri in the upper panels at an energy transfer of about 15 — 25 MeV larger 
than the experimental data. Note that these data do not include any Coulomb distortion 
corrections, which may shift the experimental result. After including the nonrelativistic 
FSI, one can see the peak position is much better represented with the error in Rl being 
negligible. 

The use of the Bonn potential in the solid curve gives a slightly different result from the 
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long-dashed curve, which uses the Franey-Love amplitudes. In this case the Bonn potential 
gives a slightly better representation of the data. 

In the middle and lower panels of Figs. 1 and 2 the transverse response is calculated. 
In this case there are two possible predictions for the same data due to the ambiguity 
about the J cc \ and J CC 2 current operators in eq. Q2.10 ). These two operators give formally 



identical longitudinal results, but differ in the transverse channel. For the nonrelativistic 
case this difference is very small. Since the A resonance is not included in this calculation the 
comparisons with the data can only be made qualitatively. From the tail of the A resonance 
in the curves, it appears that the nonrelativistic FSI predictions will underestimate the data. 

The relativistic FSI results for \q\ = 410 and 550 MeV/c are shown in Figs. 3 and 4, 
respectively. Here the predictions of Rl are very close, where the Bonn potential gives a 
slightly larger Rl than the Franey-Love amplitudes. The peak positions are accurate and 
the overall peak is well-represented. For the transverse case the cc2 or Dirac current results 
shown in the middle panels greatly underestimate the data, while the ccl current results are 
very close to the data, although the inclusion of the A-resonance degrees of freedom may 
easily alter this agreement. 

From the Rl comparisons the relativistic dynamical effects are very important in provid- 
ing good theoretical predictions of the data at these intermediate momentum transfers. For 
greater momentum transfers one would expect that relativistic effects to be an important 
and necessary ingredient for any accurate theoretical description of the data. 

The quasielastic 56 Fe(e, e') results for a momentum transfer of |g| = 1.14 GeV/c are 
shown in Figs. 5 and 6 using nonrelativistic and relativistic FSI, respectively. Here the 
PWA result for Rl places the peak positions at an energy transfer which is about 55 MeV 
greater than the peak in the data. With nonrelativistic FSI, the peak position is closer to 
the data but still about 30 — 45 MeV too high. With relativistic FSI in Fig. 6 the Dirac 
dynamical effects move the predicted peak further away to be about 40 — 50 MeV higher in 
to than the data. In this case the shapes of the peak in Rl are well- represented, although 
a bit too wide, while the peak position is not as accurately placed. The relativistic FSI 

8 



calculation gives a smaller amplitude peak than the nonrelativistic FSI calculation and is 
closer to the circled data, which take into account in an approximate fashion the Coulomb 
distortions. It is interesting to note that the relatively large errors are not truly able to 
discern clearly between the relativistic and nonrelativistic FSI calculations. The Coulomb 
corrected data appear to favor the need for relativistic FSI. The predicted peak position in 
Fig. 6 is actually closer than the most sophisticated results of Ref . |§ , which places the peak 
position at an uj of about 85 MeV higher than the data. It is clear from Fig. 6 than FSI 
are very important effects at even this high momentum transfer, although one would like to 
obtain a more accurate prediction of the peak position. 

In some sense in Fig. 6, 40 MeV does not seem like a large number, but on the scale 
of Figs. 3 and 4 this value becomes significant. The incorrect peak position seen in Fig. 6 
cannot be interpreted in terms of an average binding energy shift, since this would also 
affect the intermediate energy range results by the same amount, where the peak positions 
are accurately reproduced. The shift in the peak position must arise from a dynamical 
effect. There are a number of possible candidates for such an effect, such as the restoration 
of current conservation, meson exchange effects, ambiguities in the off-shell structure of the 
electromagnetic current, FSI effects not included in an optical model description, better 
treatment of the Coulomb distortions or even the need for a fully causal or Lorentz invariant 
description. With the advent CEBAF, comparisons with other more accurate experimental 
results, especially in the (e, e'p) case, may prove to be very enlightenning. 

IV. CONCLUSION 

The longitudinal and transverse response functions for the inclusive quasielastic elec- 
tronuclear scattering reaction from 56 Fe are calculated using a relativistic finite nucleus 
model with FSI and Dirac dynamical degrees of freedom. FSI are included using the optical 
model Green's function doorway formalism, which for the case of a one-body electromag- 
netic current operator provides a consistent description of the final states resulting from 
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the reactive content implied by the imaginary part of the optical potential. The optical 
model Green's function is calculated in a fully-off-shell momentum space calculation us- 
ing optical potentials derived from first-order multiple scattering theory using two different 
nucleon-nucleon interactions. 

It is found that at intermediate momentum transfers of 410 and 550 MeV/c, the relativis- 
ts FSI calculation gives a very good theoretical description of the data, reproducing well the 
position, shape and amplitude of the quasielastic peak. At \q\ = 1.14 GeV/c comparisons 
with recent data find that the shape and amplitude of the peak are well-produced, but that 
the peak is placed at an energy transfer of about 40 MeV higher than the experimental 
result. The source of this discrepancy is not clear and may provide a motivation for future 
investigations. 
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FIGURES 

FIG. 1. The inclusive quasielastic separated response functions are shown for scattering from 
56 Fe at \q\ = 410 MeV/c. The longitudinal response is shown in the upper panel. The transverse 
response function calculated using the cc2 (Dirac) and the ccl electromagnetic nucleon current 
operators is shown in the middle and lower panels, respectively. The short-dashed curves rep- 
resent the PWA calculation. The solid and long-dashed curves correspond to the nonrelativistic 
FSI calculation using optical potentials calculated with the full Bonn potential and Franey-Love 
amplitudes, respectively. The data do not include any Coulomb distortion corrections and are from 
Ref. 0. 

FIG. 2. The same as Fig. 1, except at \q\ = 550 MeV/c. 

FIG. 3. The same as Fig. 1, except the FSI are calculated using the relativistic model including 
Dirac dynamical effects. 

FIG. 4. The same as Fig. 1, except at \q\ = 550 MeV/c and the FSI are calculated using the 
relativistic model including Dirac dynamical effects. 

FIG. 5. The same as Fig. 1, except at \q\ = 1.14 GeV/c and the data are from Ref. The 
data represented by the circles accounts for Coulomb distortion effects by using an effective q, while 
the squares omit this correction. 

FIG. 6. The same as Fig. 1, except at \q\ = 1.14 GeV/c and the FSI are calculated using 
the relativistic model including Dirac dynamical effects and the data are from Ref. Q. The data 
represented by the circles accounts for Coulomb distortion effects by using an effective q, while the 
squares omit this correction. 
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Fe, for q = 550 MeV/c. 
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